AD-A115  307  FOAEXGN  TECHNOLOGY  DIV  WR I 6HT-PATTERS0N  AFB  OH  F/6  20/4 

THE  NUMERICAL  SOLUTION  OF  MACH-2EMNDER  INTERFER06RAM  ANALYSIS  0--ETC(Ul 
APR  82  P  0IN6»  S  PAN 

unclassified  FTO-ID(RS)T-OO16-02  NL 


FTD-ID(RS)T-0016-82 


FOREIGN  TECHNOLOGY  DIVISION 


THE  NUMERICAL  SOLUTION  OF  MACH-ZEHNDER  INTERFEROGRAM 
ANALYSIS  OF  SUPERSONIC  AIRFLOW  ABOUT  CYLINDRICAL  AND 
SYMMETRICAL  PROJECTILES 


by 

Ding  Peizhu  and  Pan  Shoufu 


Approved  for  public  release; 
distribution  unlimited. 


Accession  For 

NTIS  GRAAI 
DTIC  TAB 
Unannounced 
Justification- 


□ 

□ 


By - 

I  D1 St  ribut 1 on/ 


Availability  Codes 
Avail  and/or 
Special 


Dist 


FJ0-ID(RS)T-OO16-82 


EDITED  TRANSLATION 


FTD-ID ( RS )T-0016-8  2 


28  April  1982 


MICROFICHE  NR:  FTD-82-C-000524L 

THE  NUMERICAL  SOLUTION  OF  MACH-ZEHNDER  INTERFEROGRAM 
ANALYSIS  OF  SUPERSONIC  AIRFLOW  ABOUT  CYLINDRICAL  AND 
SYMMETRICAL  PROJECTILES 

By:  y^)ing  Peizhu  andean  Shoufu 

English  pages:  16 

Source:  Acta  Scientiarum  NaCuralium  Universitatis 
Julinensis,  Nr.  3,  1981,  pp.  60-71 


Country  of  origin:  China 
Translated  by:  LEO  KANNER  ASSOCIATES 
F33657-81-D-0264 
Requester:  FTD/TQTA 

Approved  for  public  release;  distribution  unlimited. 


THIS  TRANSLATION  IS  A  RENDITION  OP  THE  ORIGI¬ 
NAL  POREICN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 

PflEPAREO  BY: 

ADVOCATED  OR  IMPLIED  ARE  THOSE  OP  THE  SOURCE 

AND  DO  NOT  NECESSARILY  REPLECT  THE  POSITION 

TRANSLATION  DIVISION 

OR  OPINION  OP  THE  POREICN  TECHNOLOGY  01- 

POREICN  TECHNOLOGY  DIVISION 

VISION. 

WP.AFB.  OHIO. 

FTD  -ID(RS)T-0016-82 


Date  28  Apr  19  82 


GRAPHICS  DISCLAIMER 

All  figures »  graphics.  Cables,  equations,  etc.  merged 
into  this  translation  were  extracted  from  the  best 
quality  copy  available. 


I 


THE  NUMERICAL  SOLUTION  OF  MACH-ZEHNDER  INTERFEROGRAM  ANALYSIS  OF  SUPERSONIC 
AIRFLOW  ABOUT  CYLINDRICAL  AND  SYMMETRICAL  PROJECTILES 


Ding  Peizhu  and  Pan  Shoufu 
SUMMARY : 


In  this  study,  four  physical  conclusions  and  the  numerical  solutions 
described  in  [1]  concerning  the  Abel  transformation  in  Plasma  Spectroscopy 
involving  shock-waves  extended  to  the  Mach-Zehnder  interferogram  analysis 
in  supersonic  airflow  about  cylindrical  symmetrical  projectiles.  Regarding 
the  relation  between  the  interferometric  fringe  shift  6(x,2)  and  the 
concentration  variation  p(r,z)-  obtained  at  each  point  of  the  supersonic 
airflow  field,  we  obtained  three  physical  conclusions  similar  to  those 
described  in  [1],  In  addition  we  provided  the  interferometric  fringe  shift 
5(x,i)  obtained  from  observations  and  the  coefficient  formulas  of  the  entire 
transformation  necessary  to  calculate  the  concentration  variation  p(r,z)  -p^ 
at  every  point  of  the  airflow  field. 

We  borrowed  from  source  [2]  the  airflow  field  density  computation 
method  provided  for  the  Mach-Zehnder  interferogram  applied  to  supersonic 
symmetrical  projectiles,  and  obtained  satisfactory  results.  Also,  a 
simplified  method  for  computation  is  given  for  the  case  of  additional 
shock-waves  present  in  the  airflow  field. 

I .  INTRODUCTION 

The  question  of  the  airflow  about  supersonic  cylindrical  symmetrical 
projectiles  is  old  and  of  sizeable  importance.  The  Mach-Zehnder  interferogram 
analysis  (*)  has  been  for  many  years  an  important  and  practical  research 
tool.  Since  the  invention  of  digital  computers,  extensive  applications 
occurred  in  the  M-Z  interferogram  analysis.  The  method  of  interferogram 
analysis  consists  in  measurements  of  the  interferometric  fringe  shift 
6(x,z),  the  research  of  the  concentration  variations  p(r,z)-p^  at  every 
point  of  the  airflow,  using  the  Abel  transformation;  thereby  we  obtained 
the  curve  of  density  in  plasma  spectroscopy, 

(*)  below:  *’M-Z  interferogram”  or  ’’interferogram”. 


In  the  past  scientists  have  applied  the  integral  of  Stieljes  to 
obtain  numerical  solutions  of  the  above  problem  [3],  But  experimenters 
are  still  not  familiar  with  the  application  of  the  numerical  solutions  of 
the  integral  of  Stieljes;  further,  the  results  computed  thereby  present  some 
ambiguities  in  this  case  of  shock-waves.  In  [1],  we  discussed  some 
characteristics  relative  to  the  usual  simplified  solution  of  the  Abel 
equation.  By.  using  the  Abel  equation  with  a  radiation  frequency  I(x)  in 
plasma  spectroscopy  involving  shock-waves  and  an  emission  coefficient  (E(r), 
we  demonstrated  that  E(r)  can  be  computed  from  the  usual  Abel  transformation, 
and  we  obtained  the  formula  E(r)  for  the  computation  of  the  amplitude  of 
the  shock-waves. 

This  is  why  we  suggested  in  [1]  a  new  computation  of  the  Abel 
transformation  in  plasma  spectroscopy  involving  shock-waves. 

In  the  M-Z  interferogram  analysis,  equations  noted  as  6(x)  and  p(r)-p^, 
the  radiation  frequency  in  plasma  spectroscopy  I(x),  and  the  emission 
coefficient  E(r)  are  mathematically  identical.  Therefore  it  is  possible  to 
extend  directly  the  conclusions  and  methods  of  computation  obtained  in  [1], 
to  the  Abel  transformation  of  interferogram  analysis  about  cylindrical  and 
symmetrical  supersonic  projectiles.  The  method  of  computation  is  simple. 

The  physical  theory  applied  to  the  shock-waves  is  equally  understandable. 

In  the  second  part  of  this  study,  we  describe  the  Abel  transformation 
of  interferogram  analysis  about  cylindrical,  symmetrical  supersonic 
projectiles.  In  addition,  we  extended  the  conclusions  obtained  in  [1] 
concerning  the  Abel  transformation  in  plasma  spectroscopy  involving 
shock-waves  to  the  present  case.  In  the  third  part  we  have  listed  numerical 
solutions.  In  the  fourth  part  we  discuss  the  case  of  an  additional 
shock-wave.  The  fifth  part  is  a  possible  plan  for  a  simplified  computation 
in  the  case  of  an  additional  shock-wave. 

II.  THE  ABEL  TRANSFORMATION 

Figure  1  illustrates  a  point-headed  supersonic  projectile  in  flight. 
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We  suppose  that  conditions  are  such  as:  it  is  cylindrical  and  syirmetrical, 
and  that  there  are  one  or  several  tapering  shock-waves  iroving  towards  the 
back.  Section  1-1  is  illustrated  on  Fig.  2:  the  hatched  circle  in  the 
center  represents  the  cross  section  of  the  projectile;  the  circumference  c 
represents  the  cross  section  of  the  shock-wave.  Outside  c,  the  airflow  has 
not  yet  been  affected  by  the  wave's  interference.  The  density  and  the 
index  of  refraction  are  respectively  and  n^;  within  c,  since  we  are 
dealing  with  a  cylindrical,  symmetrical  projectile,  the  density  and  the 
index  of  refraction  only  reach  the  distance  r  related  to  the  center  0,  and 
are  respectively: 

P=P  (r) ,  n*n(r) 


Fig.  1.  Illustration  of  the  cylindrical,  Fig.  2.  Illustration  of  photoline 
symmetrical  supersonic  projectile  according  to  cross-section 

in  flight.  1-1. 


Fig,  2  illustrates  the  interferometric  fringe  shift  corresponding  to 
section  1-1.  The  value  6(x)  at  point  x  is  the  light  crossing  section  1-1 
refracting  the  superinposed  effect  of  one  beam  going  from  -y  to  y^  and  away 
from  the  y  axis  for  each  value  of  x. 
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According  to  the  usual  hypothesis,  the  interferometric  fringe  shift 
6(x)  and  the  refraction  index  difference  n(r)-n^  or  the  density  difference 
oCr)-o^  satisfies  the  equation; 

<3(x)  5=-^L|  *  (Fi(r)-n,)dy=  \ 

in  which  K  is  the  Gladstone-Dale  constant,  is  the  wave  length  in 

monochromatic  light  and  total  vacuum  necessary  for  obtaining  the  interferogram. 

2  2  2 

Focusing  on  the  cylindrical  characteristics  and  the  value  of  x  +y  =r  ,  the 
above  equation  becomes; 


2K  r  (f*(r)  -  /I„)r  . 

dr 


(1) 


We  assume  that  the  shock-wave  is  a  simplification  of  the  first 
categor>"  of  o(r)-o^  (or  n(r)-n^).  If  c  is  the  most  external  shock-wave,  then 
outside  c,  o(r)-D^=0.  And  therefore  Eq.  (1)  can  be  written; 


Six) 


2K  VHp(r)  - 

A-  J// 


(2) 


in  which  R  is  an  arbitrary  constant  superior  to  c.  Besides  the  constant 
2K/X*,  Eq,  (2)  and  the  Abel  equation  encountered  in  the  study  of  plasma 
spectroscopy  involving  shock-waves  are  similar.  Bennett  and  his  colleagues 
have  studied  Eq.  (2)  in  [2];  in  [1]  we  discussed  some  characteristics  of  the 
solution  of  the  Abel  equation.  If  we  apply  them,  to  Eq.  (2),  we  obtain: 

1)  on  the  M-Z  interferogram* s  arbitrary  cross  section  1-1,  the  fringe 
shift  is  continuous. 

2)  If,  on  section  1-1,  the  shock-wave  c  belongs  to  the  first  category 
interrupted  p(r)  and  p’(rj,  then  inside  the  shock-wave  we  have: 


S  (jc) 


ZKx  f*  p'  (r) 


dr* 


ZKx 


J 


(3) 


especially  when  x-^. 


2)  If  the  fringe  shift  is  srooth  outside  the  peaks  c,  then  d'W 
exists;  when  x-^c,  d'(x)^  ^/J  c-x  t  then  shock-waves 

appear  outside  the  airflow  circle  r*c,  on  the  cross  section  1-1; 
the  density  saltatory  value  of  the  shock-waves  will  be 

=  (5) 

KJ  2c 

further,  on  [c^,c]  and  [c,R] 


Tlaese  conclusions  are  also  applicable  in  the  case  of  a  single 
shock-wave. 

III.  NUMERICAL  SOLUTIONS 

We  will  first  discuss  the  case  of  c  with  a  single  shock-wave.  For  a 
convenient  application  to  the  case  of  several  shock-waves,  ’’The  p(r)-p^=0 
external  to  c”  is  transformed  into:  ”on  [c,R],p(r)  and  p’Cr)  are  continuous 
when  r^R,  p(r)-p^=0.”  Evidently,  the  former  is  a  particular  of  the  latter. 

With  [0,R]N,  we  have  the  following  partition; 

o  *  X,  .  <  ..,<0  =  X,  <x,^ ,  <  =  R 

or  o  =  *.<- •<c,  =*•<*•♦,<■  •<*.=«<****<"<*''“^* 

We  also  note:  A  =  x^^i-x,,3g  =d(Xg)  ,  According  to  conclusion  2 

obtained  in  section  II,  the  derivative  6’(x)  of  the  interferometric  fringe 
shift  tends  toward  the  moment  of  the  shock-wave  with  the  value  1/^c^x 
and  toward  the  infinite  fro®  inside  the  shock-wave.  Consequently,  one  more 
term  is  used  in  the  expression  of  the  shock-wave,  1  ic*-x*  ^  for  the 
breakdown  of  the  couqparable  6(x).  On  extract 

^#-#(x)  «X4fx  +  ^x*  +  /j"P^rjT" 

the  comparable,  5(x),  zero. 
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IVe  can  state 


3,y-  2,y-  1 


including  a,  b,  c. 


/  = 


3^11-  I  ”  33,  .f 


(8) 


Tlien  extracting  from  D^3j  2v-l  ~3J  4v-V  +J  6v-9  . 
C*<,*»  +  iD(/  =  w+ l.w  +  2,  — ,w-3) 

P,  (*)  =a,  +6,jf +  C(Jf*  +d,x*  +//‘c^- jc*‘ , 

IVhen  Cq^O,  we  extract  from 

/*•  (X)  =  a,  +  b^x  +  c^x*  +  IJ  c*-x*  , 

IvTien  c^=0,  we  extract  from  [O.Xj^] 

Ptix)  =a,^.b,x  +  ctx*  ^  d,x»  +  IJ-pzr^, 


The  comparable  interferogram  fringe  shift  <S  (x) ;  in  every  small  interval 
external  to  the  shock-wave  +  =t;,y  +  ]  ,-.,7^- i)  we  use  the 

ternary  or  binary  polynomial: 

Pfix)  =a,  +fr,x  +  c,jc»+c/,jc*, 

Pfix)  =a,  +6#jc  +  c,x»  (;  =w,iV-]) 

as  the  comparable  interferometric  fringe  shift  6Cx).  See  [4],  By  using 
again  the  value  5^,  we  can  state  the  various  coefficients  a^,b^,c^.  According 
to  the  conclusions  obtained  in  section  3 


p(r,)-p.=  -  [*  dx^--^y2r*‘  --^L^dx 

m  7  f 


By  using  P'^(x)  as  stated  above  to  substitute  6'(x),  we  get: 

A*  P'kW_dx 

P  (r , )  -  P«  =  -  i:  j^rTi*  * 

After  computation  of  each  integral,  we  have  the  parallel  systematization: 

.V-  I 

P(r*)  =p,  +  iE  o,.*  d,, 

/T 

.V-  I 

I  ,a;4»  -  2ey,y  +  1  ,  —  eiV  «  ])  I 

/r  k  i  -  t 
.V-  I 

P<r..,)  ^P„^l  H  a,.,.,d», 

^  -  .1 


C9) 


We  will  note: 


C* 

D-  :?/?,  - 

/  //,- 

/.(» 


t'm 


yUw) 


\2^K'  -J2v~t  JK-J  IV- I  ./f,=Jgv^a  f  Hv-  IBf 

+  /?,  ,/r  =  2/?,  -  /},,/•■  r  zn,  -  n,  ,g  =  35,  -  2/?, .//  =  n,  -  zb,  +  /?„ 

;/^  +  3^,-//.,  /-■(/./<) ---((i/e— 2  +  3;-')  In  (&+  -^kj  , 


T(i  .k) 


arctuj 


y  ^  Jr'  -  j  •• 

J  (:--2)  - 


y- -(/(-(-  !)- 

(■4,  H  J  iT_7,. 

/  (;■'  -  Jv-  -  j-  -  J  (r--2)’-/», 


I  -t- 1 2rJ  l- '  -  1  -arcsin  ^  + 

*  ( J  y- - 9  -  3i  y*  -  4  +  3j  w"  -  I  -  y)  •/■'(rt.  1 )  ~ 

12 

-  (O.rJ  y“  -  f  -  2J  y»  -  I  +  )  ..iw)  •  l''2 

^  J  iP  -  1  -  2V  y- -  1  +  +  7" (0 . 1 ) + 

+  12  +  +6ju*-fc*  -  4J  yf- (/e- 1)* 

+  Jy— (/e-2)-’ )  •F(0./<)+r(0,*)[-  I  •  F  (0  ,v  -  2)  + 

+  {(y-2.y)5,  -2<t'-2)5,+  (y-  l.r>)5,}In(2w-4)-//*(y-2), 

+  (w  +  2)*  -2(a»+  1)  J  y*  -  (u;  +  ! )  ^  + 

+  (a;  +  1 ,5)7  J'*  *  w*}  Inu;  + 

+  !  ^7  y*  -  (w  ^  3)  ’  -  37  y*  -  (w  +  2)  *  +  37  y *  -  (a;  +  I )  -  - 


-  7  V*  -  w’l 


•  +  1  )  + 


y  --f " 


+  E  (A +  i)*+6J  «•-<*'  ■ 


-4Jw»-  (A-  1)*+JW*-  (S'^2)*)*7'’(».*)  +7’(w,A)}- 
-  ^/.F(w,w-2)  +{(w~2.5)^,-2(»-2)5.+  <w-1.5)flt}* 
•  In(w-2  + J  (»-2)*-w*  )  (W“2)*- w^» 
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A<;)  \  -  0,:))  J  -K  I)  »  -  2/V  t;* -7*  +  O'  +  0.5)7 7 ^  -  J )  -  }  !«;  + 

"  i2{-'y*~T;>2)  *  -  ••JJy*-  (;■  + i)*  +  3J17^*  - 

-  J  y-  -  (;  -  I)  ‘I  /•' (;■  .;■)  +  arct^|,'  ~  +J>.i  _  -1  + 

lli!  1 

^  12  '  rj  (-  -  (*  +  2)  *  -  -1-'  v'^~7kT Ip  +  6-'  t-'  -  - 

t  »  ♦  j 

-  t7  y»-  (/;  -  j)»  +  +7’(;,/f)|- 

-  /•/'’ (/.'p- 2)  +|(y-2.5)^,-2(y-2)^,  + 

+  .ro  /?,  |ln  (y-  2  +  7  (y-  2)"*  -  -7  <y-  2)  ’  ^  « 

(j  =  w  -h  I, •••.£/"■  3). 


Each  coefficient  a.  .  in  the  fornula  of  the  transformation  (9),  is 

^  ^  -  ■  -  ^ 
expressed  below  (^supposing  w+o<v<N-o} 

^In2/  \viicii  iXr^' 

n  ' 


,125  +  ' 


j  1 2C**  I  i.iitj-*-  ^ 

(C*(F(tt;  ^a;  +  2)  *  F  (yj  ,u/ +  O  -  5(2u;  ^  when  \\>\ 

,12C*|-I.5+ 

\c*{F(w,w  +  3)  -  iF  (ii\zv  -^2)  +  3F  -t-  !)  4*  2  J  (u'  1 )  w  ^  I 


-I2C*|o.375-  ^3^^”-} 

Ir* 


when  w  =  ;i 


C*{F(w,w  +  4)  -  4F(u;,a;  +  3)  +  fiFCty.iu  +  2)  -  3F(tu,:y  +])- 

-  6(2iy  +  1) Injjy}  when  ;7^i 

=  C*{F(; ./«  +  2)  -  4F(;\&  +  I)  +  oF(j\k)  -  iF(j\k-  O  +  F(j -  2)}» 

(;  =  a;,A  =  w  +  3,u^  +  4,  —  ,^-  4)  j 
(;  =y;  +  Kw  +  2,  —  .y-  6*^  =  ;  +  2,;  +  3.  — 4) ; 
ij  =  w^w+  1  »•••  .y  +  1  =  y  +  3  .y  +  4  .•••  •^V  -  3)  I 

(;  =  y  +  2,y  +  3,*«»,:V -  5,^  = ;  +  2,;  +  3,—  ..V-  3) , 

=  C*{F(;J  4-  I)  -  F(jJ)  -  6(2;  +  ))ln;} , 

(;  =  uj  +  1  ,a>  +  —  3,y  +  1  ,y  +  2,  —  ..V  -  2)  j 

<if.#  =C*{F(;,;  +  2)  -  iF(;,;  +  1)  ^zFijJ)  +24;In;}, 

(;  =  tv  +  1  ,tv  +  2,  —  ,y -  4  .y  +  I  ,y  +  2 ,•••  ..V  -  3) : 

--  C*(FO'.;  +  3)  -  4F(;\;  +  2)  +  fiFr/,;  -v  i)  -  ZF ij  J)  -  G(2;  -  Dln;^ . 

(;  _  iv  4-  )  ,tu  +  2  ••••  -  5  .y  +  1  .y  +  2,***  •  A"  -  0  i 

tr-s  =C*{F(;,y-  5)  -  4F(;,y-  4)  +  sFO.y-  3)  -  3F(;,y-  2)  -r 
+  6(2y  -  3)  In  (y  -  2  +  7(y -  2)  *  ^  )  “  J27  (y  -  2) '  -  ^  “ 


(/  tu.tv  +  I 
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,y  -  5) 


y-2)-’-;’ 


a,..--, -C*{F(;.u-4)  -4F(;.y-3)  +  3F(;  .y  -  2)  -  21  (t>  -  2)  <  In- (0-2  + 

+  J  (y-2)  =  -;=  )+24J  }- 

-F-/ (;■)  -^Sarctu;  j  J’ 

(j  -  w.ui  +  I ,v  -  1)  i 

a,„.  ,  -C’{F(y.y-3)  -F(;.-J-  2)  +  f.(2t' -  5)  In  (i;  -  2  +  J  (i'-  2)  =  -;’)  - 

-  12J  (t>- 2)’ [[- D+:’,B,-  <'7]-L(j)  +(7-J  (j)  + 

R 

'I- 3^(7)  -  ^arctt’^  '  ,  (;  u;  .’c» -V  1  ; 

J 

a, F(;\y  J-  2)  -  F  f;  .i'  +  I )  -  r,(2z^  “  y  '■  ^  “ 


-  '  -|<F-  />’,  -  //)./.(;)  »-  //-y  (y>  -  .J(»  +  arcli;  j  („_!>*  j  _  ya  .  ’ 

(;  =  a;,Lt;+  I ,y-  3)  ? 

(i,  ,-C*(-2F(y-4.y-n+r.Ffy-i.y-3)-:)F(i'-l.y-2)- 

-  li  (2:’  -  iO  !>'  <i’  -  t)  +  ('i(2:’  -  .’I) !  >'  (i'  -■  2  +  2J  1  ~  2  I J  -  :i  )  “ 

-  I)  +AVJ-  D-arcfi;  i 

D  \v  2  '  2Jv-i 

^  C*{3Fd’-  3.ti-  3)  -  ?.F  (V  -  2 ,1’  ~  2)  +  0(2^  -  3)  In  (i- -  2  +  J  2^  -  5  >  - 
-  I2j2t'-r)  +24{t'-3)ln(:,-3)}-  '  {(l  “  F .  • -  3)  + 

+  F-y  (u- 3)  +  -  3)  -  arct?  }( 

J  2t'  -  .1 

F*{3F(y-3,-J-2)  -3F(n-3.y-3)-6(2u-7)ln(:-~3)  - 

-  2  I  (y  -  2)  In  iv-Z-^J  .')  )  +  2  I J  2u  -  a  I  “ 


'  '  I  ( -  3  +  ^ )  F . .  f  -  3 )  -  F  •  y  ( .V  -  .! )  - 


B  ) 

-  3^(y- 3)  +  3arct(;  - — Vi 
J2V-5  I 


a  -  '2C*f  .T  F,.Fx  5  + 

^  \  2  V*  2  '  f-2 

a.. . ,-L2f--{-|.,.fl.«£-2D,t[D-3B,*  f  <3£- zD) ]l.  } 

......  - 

+  C*{F(t;-2,a  +  2)  -F(w-2,w+  l)-6(2«  +  3)ln(t;  +  fl,)}, 

a . 
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0.-...-.=  •|,--{--|7«  +  5.(3£--D)+  [3-|j-  (3£-D)j5..1n-^^}, 
a, -  fi .  (3E  -  2D)  +  [ D -  3fl .  +  (3E  -  2D)  ]ln  , 

a.- ...  =  T  (E-D)  ]ln  "  + }  + 

+  C*{F(y-  1,0  +  2)  -/^(w-  1.W+  l)  -6(2o  +  3)ln(y+5,)}, 

«...  +  .  =C*{F(;,o  +  3)  -  4/!'(;,v  +  2) +3F(y,o+ 1) -2W*'^»'  + 

+  24(0+  l)ln(o  +  ^  0*-;*)  .  O  =  «>,«>+  l.—  .y)  i 
.  ="C*{F(;,o  +  4)  -  4F(;,o  +  3)  +  6F(;,o  +  2)  -  3F (; .o  +  ])  +  12jJ’-7»- 
-  6(2o  +  1)ln(o  + J  [I*  - 7») } ,  (;  =  w,UF+  1  .•••,o). 

‘J.-.-.  =0i 

a,..,  =  C*{  -  F{y,o  +  1)  +  F(y,o  +  2)  -  G(2o+  3)  Ino} . 

=C*{F(;,.V-  4)  -  4F(>,yV  -  3)  +6F(;,A(-2)  -  3F(;,/V  -  i)  + 

+  6(2.'V  -  l)ln(iV  +  .7  ;V*  -  >*  )  -  I2j  (;  ^  w,io+  l,  —  ,iV  -  4) , 

1  C*{F(;./V  -.3)  -  \E(i,N -2)  +.3F(;.7V-  I)  - 

-21(.v-  Dliu.V  +J  ,V*-F  )  +24.7  iV'-j’  }  (;■=-«;,«'+  — ,^V-  .3) , 

C*{-3F(.'V-3,7V-3)  +6F(/V-3,iV-2)  -3F(.V  -  3,7V-  l)  - 
-6(27V-7)ln(7V-3)  +6(27V-  l)ln(.V+.7  67V-9>-  12.7  67V-9  }. 
iv.,.v-*  =C’*{3F(7V-2,iV-2)-3F(7V-2,7V-  1)  +  6(27V  -  1)  ln(iV  +  2j  7^^) 
-24jiV^l  +24(iV-2)In(:V-2)}; 

“v-j-x-i  =C*{3F(.V-2,7V-l)-3F(7V-2,A7-2)-21(iV-  !)• 

•  In  (TV  +  2777"-  1)  +  48  777^ -  6(2iV  -  5)  In  (TV  -  2) } , 

a,v.i.,v-5  -C'*|6(2.V  -  _  12  727V-  1  j, 

"  C*{  -  24(7V  -  1)  In  ^TV  -  i  +2  4  727^-l}. 

According  to  conclusion  2  obtained  in  section  II»  and  Eq.  (7),  it  is 
possible  to  establish  the  approximate  rebound  equality  for  the  density  of 
the  shock-wave. 
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p(c.)  -  p(Ct)  =  X*ll2K, 


(10) 


in  which  1  is  computed  as  in  (8).  Especially  when  r^c,  and  in  the  case 
when  p(r)-D^=0 


P(c.)-P,  =  A*l/2Ar.  (11) 

Supposing  the  measurements  of  the  interferometric  fringe  shift  give 
the  following  readings: 

/(X) 

<5(x)=( 

^f(x) 

and  if  in  addition  f’(x)  is  continuous  with  and  f ’ * (x)  a  fractionate 

succession,  we  have  proven  (in  a  further  study  on  the  computation  of  errors) 
by  the  above  calculation  of  the  density  of  the  airflow  p(r^),  the  existence 
of  an  error  such  as: 


ilP(ri) 


(12) 


in  which  the  constant  M  and  c,  R  and  are  related,  and 
the  positive  number  a  equivalent  to  the  wave  h,  being  extremely  small,  can 
be  arbitrarily  small. 


I  -  p(c^)) 


max 


\r  Ui)  -/" (A'j)  |// ,  (13) 


The  constant  M*  and  c  are  related.  From  computation  in  (13)  it  is 
obvious  that  the  approximate  calculation  (10)  is  rather  precise. 

From  Table  9  established  in  [2]  regarding  the  interferograms  of  spherical 
projectiles,  we  have  estimated  the  values  of  6(x,2),  the  interferometric 
fringe  shift,  and  used  them  in  the  computation  with  the  method  described 
above.  Fig.  3  shows  the  isograph  of  the  comparative  density  values  o/o^ 
inside  the  shock-wave.  Fig.  4  is  the  reproduction  of  Fig.  10  borrowed  from 
[2].  Table  1  gives  the  concentration  saltatory  values  of  the  shock-wave 
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Table  1. 


I 

I 

i  1 

i  ff  i 

i  V 

0.61 

1 

0.26 

0.19 

0.14 

1.60 

!  1.60 

1.44 

1.38 

1.33 

IV.  CASE  OF  SEVERAL  SHOCK- WAVES 

When  a  point  headed  cylindrical  sxxpersonic  projectile  is  in  flight, 
two  or  more  shock-waves  can  be  observed.  As  illustrated  in  Fig.  1,  cross 
sections  2-2  and  3-3  correspond  respectively  to  two  or  three  shock-waves. 


With  two  shock-waves  as  an  example,  supposing  that  the  wave  originates 
outside  r-c^  and  r=C2,  the  projectiles'  cross-sectional  radius  is  c^: 

^0  'Cc, 


we  extract  the  isometric  partition; 


'  c  ^ 


In  the  expression  w.w  +  J  -  1)  we  extract  the 

binary  or  ternary  polynomial  comparable  6(x)  with  the  additional  term 

l,Jc;-x’  •  expression  C** .*>♦  i3  o’ =  »Wi  +  1  *•••»''. r  1)  we  extract  the 

binary  or  ternary  polynomial  comparable  6 (x)  with  additional  term  l,J  c|- x*  • 
By  repeating  the  process  described  in  section  III,  we  can  definitely 
extract  1^^,  I2  and  the  coefficients  of  each  polynomial,  introducing  them  in 


P(r()  -p,  = 


UK 


d.X  , 


(/  -  u;,w+  1) 


which  gives  us  the  parallel  systematization: 
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V  «  -  I 

/J(r».)=p,  + 

*t-* 

/>(ry)=iP,+  ^a(.»<5t,  (/■  -  iu+ l,  —  ,y,  -  2.!»,  +  - 2) ; 

-  t 

*•*  I 

P(r,  -P,+  a,,_ 

»-*, 

’f* 

P(r.  )-  P»+  2  a(.t<5fc,  (I  y ,  -  1  .y.  -  1) . 

k  - 1  -  » 

in  which  each  equality  of  transformation  coefficient  a^  is  comparable  to 
the  derivation  of  section  III, 

On  the  external  shock-wave  c^,  the  concentration  saltatory  value  is: 

p(c,_)  =.  ^ 


On  the  internal  side  of  wave  c^,  the  value  is: 

-  Pd  +  A  In/ jA  • 


in  which: 

-  Q  5  Adp  ^  ^ .  +  (\,  ,  <  j  , 

On  the  internal  side  of  shock-wave  from  the  concentration  saltatory 
values 

~p(c,^)  aU,/2K 


and 


■'  9  '  ' 


we  obtain  the  concentration  rebound  values  of  the  internal  side  of  the 
shock-wave: 


n(c ,  ) 


2K 


in  which: 


/|  -  j "  *  *"^1* ,  -  *)  • 

=  3j  2y,  -  I  -  3J  4y,  ITj  r  J  -  9  . 
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V,  PLAN  FOR  A  SUIPLIFIED  COMPUTATION 


\ 


We  take  as  an  example  the  case  of  two  shock-waves.  The  partition  is: 

Previously  we  only  considered  one  shock -wave  c^.  As  in  Eq.  (9)  in  section 
III,  we  extr  v=v^^,  w=Vj^,  to  compute  the  value  of  the  concentration 

1 ) »  •••,  P(r,.  , ) 

IVe  distinguish  the  values  of  concentration  of  the  shock-waves  c^^  and  c 


The  rebound  value  of  the  airflow  concentration  of  the  shock-wave  c^  is 

P(c,-  )-Po  =  A*/5/2/C, 

in  which  1^  is  similar  to  section  IV.  Secondly,  we  will  review  the  case  of 
one  shock-wave  Again  according  to  Eq.  (9),  we  extract  v=v^,  N^v^,  to 

compute  the  concentration  value: 

P(r»),*—,  P<r,.,),  p(r,) 


We  distinguish  the  concentration  values  in  the  shock-wave  z^  and 
outside  of  it; 

P(rw),  •••,  p(r,^.,),  p(c,^ 

The  airflow  concentration  values  are: 

P<c,.) -P(c.*)  =<1*/,/2A:, 

in  which  1  is  similar  to  the  value  obtained  in  section  IV. 


The  method  described  above  can  repeatedly  be  applied  to  the  case  of 
several  shock-waves. 
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Abstract 

The  method  of  numerical  solution  on  four  physical  conclusions  and 
“The  Abel  Transformation  of  Plasma  Spectroscopy  Involved  Shock  Wave** 
obtained  by  the  same  authors  in  ^11  is  extended  to  Mach-Zehndcr 
interferogram  analysis  of  supersonic  airflow  about  cylindrical  symmetric 
projectiles.  Three  similar  physical  conclusions  about  the  relation  between 
the  interferometric  fringe  shift,  measured  from  M^Z  inferferogram 

and  the  density  difference,  p(r,^)*p„,  at  each  point  in  the  supersonic 
airflow  feild  are  obtained,  and  all  transformation  coefficient  formulas  used 
to  calculate  the  density  difference,  p(r,c)  -  according  to  the  interfe¬ 
rometric  fringe  shift,  (5(-v,r),  have  been  given.  As  an  example,  by  use  cif 
experimental  data  given  in  C2^,  airflow  field  of  cylindrical  symmetric 
projectiles  is  calculated  and  the  better  results  have  been  obtained,  A 
simplied  calculation  method  which  deals  with  many  shock  waves  existiui; 
in  airflow  frild  has  been  given. 
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